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Abstract. The adiabatic index Fi for perturbations of dense matter is studied under various physical conditions 
which can prevail in neutron star cores. The dependence of Pi on the composition of matter (in particular, on the 
presence of hyperons), on the stellar pulsation amplitude, and on the baryon superfluidity is analyzed. Timescales 
of damping of stellar pulsations are estimated at different compositions, temperatures, and pulsation amplitudes. 
Damping of pulsations by bulk viscosity in the neutron-star cores can prevent the stars to pulsate with relative 
amplitudes J> (1 — 15)% (depending on the composition of matter). 
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1. Introduction 

The adiabatic index Ti = (rib /P)(dP/ drib) for density 
perturbations determines the changes of pressure P asso- 
ciated with variations of the local baryon density rib (e.g., 
Shapiro & Teukolsky 1983). It enters the equations gov- 
erning small-amplitude neutron-star pulsations (Thorne 
& Campolattaro 1967, Thorne 1968, 1969) as well as the 
criteria of stability of cold relativistic stars (Meltzer & 
Thorne 1966, Chanmugan & Gabriel 1971, Chanmugan 
1977, Gourgoulhon et al. 1995). For the time-dependent 
perturbations, like neutron-star pulsations, Ti has to be 
calculated taking into account the slowness of various 
equilibration channels in dense matter. The same factors 
strongly affect viscous damping of neutron-star pulsations. 

In the present paper, we study three main factors 
which regulate T\ and the viscous damping of pulsations: 
composition of dense matter, pulsation amplitude, and su- 
perfluidity of baryons. Equilibration processes for various 
compositions of neutron-star cores are studied in Sect. 2. 
In Sect. 3, we calculate the adiabatic index under various 
conditions of density, composition, and temperature. We 
separately study two different regimes: first, the regime 
with the perturbations of chemical potentials of particles 
much smaller than T, and second, the regime with the per- 
turbations much larger than T (we use the units in which 
the Boltzmann constant fee = 1); we also consider the 
effect of baryon superfluidity. In Sect. 4 we estimate the 
timescales of the viscous damping of density pulsations 



under various physical conditions which can be realized in 
the neutron-star cores. Finally, in Sect. 5 we summarize 
our results and briefly discuss the problems which remain 
to be solved. 



2. Relaxation processes in neutron-star cores 

It is widely believed that dense matter of the outer 
neutron-star core (0.5 no rib & 2 no, where no ~ 0.16 
fm -3 is the number density of saturated nuclear matter) 
contains neutrons (n), protons (p), electrons (e), and - if 
the electron Fermi energy exceeds the muon rest energy 
(that occurs at % ^ no) - also muons (/i). At higher 
densities in the inner core, hyperons may appear, first of 
all, S~ and A hyperons. All constituents of matter are 
strongly degenerate. We restrict ourselves to this compo- 
sition and will not consider exotic models of matter con- 
taining pion or kaon condensates or free quarks. We will 
study a neutrino-transparent neutron-star core formed in 
about 30 s after the neutron star birth in a supernova 
explosion. We ignore thus the early stage of the neutrino 
opaque core in a protoneutron star (the adiabatic index 
in such matter was studied by Gondek et al. 1997). 
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2.1. Strong and weak-interaction processes 

Full thermodynamic equilibrium in dense matter is estab- 
lished only after partial equilibria are achieved in all equi- 
libration channels. In our case, there may be five types 
of channels [labelled as (a), (a'), (b), (c), and (d)], with 
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drastically different relaxation times. These channels, de- 
scribed partly, e.g., by Yakovlev et al. (2001) and Haensel 
et al. (2002), are as follows: 
(a) Strong and Coulombic elastic collisions, 



These estimates are valid if only the pulsation amplitudes 



BB' # BB', £B # £B, ££' ^ ££', 



(1) 



where B and B' are baryons (nucleons B=N or hyperons 
B=H), while £ and £' are leptons in dense matter (elec- 
trons or muons). These processes do not change particle 
fractions. 

(a') Strong-interaction processes conserving strangeness 
but changing baryon fractions, 



nA ^ pE~. 

(b) Modified Urea processes like 
nN->N'p^, pN£->nNV £ , 



(2) 



(3) 



where vi and i>i are neutrino and antineutrino. Similar 
processes can also proceed with one or both nucleons re- 
placed by hyperons. These processes can change the frac- 
tions of nucleons, hyperons, and leptons. 
(c) Direct Urea processes involving nucleons and hyper- 
ons, 



n — > p£vt, p£ — » m>i; 
£ _ — > n£&e, n£ — > Y^v^; 
A — > p£i>i, p£ — ► kvt, 



(4) 
(5) 
(6) 



which, similarly to processes (b), can change the fractions 
of nucleons, hyperons, and leptons. 

(d) Nonleptonic collisions of baryons changing 
strangeness, e.g., 



NA ^ Nn, nn ^ pS" 



(7) 



they can change the fractions of baryons, but do not affect 
the fractions of leptons. 

2.2. Relaxation times 

Relaxation times in channels (a) and (a') are about 
(10- 16 -10- 19 )/T 9 2 sor shorter, where T 9 = T/10 9 K (this 
can b e deduced, e.g., from the results of Flowers & Itoh 
1979| ). Therefore, these channels guarantee instantaneous 



(on hydrodynamical time scales) partial thermodynamic 
equilibration to the Fermi-Dirac distributions of particles 
of all species j with a temperature T and chemical poten- 
tials Processes (a') allow for additional equilibration 
of A, S~, and nucleons, which leads to the equality of the 
chemical potentials /is + fi p — fi n + ha- The relaxation in 
channels (b), (c), and (d) (with relaxation times tm, td, 
and th) is associated with weak interactions and lasts, 
therefore, much longer; it will be important for our analy- 
sis. In a nonsuperfluid matter close to full equilibrium the 
relaxation times can be estimated as (e.g., Yakovlev et al. 
2001], Haensel et al. |2002[) 



2 month 



20 s 
1 d 



1 ms 

1 9 



(8) 



djJLj = — fij | are very low, Sfij <C T, where /i^' 1 is the 
fully equilibrium value of fij. We will call this regime as 
the regime of subthermal pulsations. 

According to Eqs. (@), the slowest are modified Urea 
processes which involve neutrinos and largest number of 
strongly degenerate particles. The most rapid are non- 
leptonic processes changing strangeness; they go via weak 
interaction but involve smallest number of degenerate par- 
ticles and no neutrinos. While estimating th and viscous 
dissipation parameters in hyperonic matter from the re- 
sults of Haensel et al. ( [2002 ), we set the phenomenological 
hyperon reaction parameter \ = 0.1. All the parameters 
in hyperonic matter are rather uncertain and can differ 
from the presented values within a factor of ten or even 
larger. 



(0) 



3. Adiabatic index 

Clearly, one can divide the neutron star core into three 
zones with different relaxation properties. The outer zone 
(zone "M") extends from the core-crust interface, rtb ~ 
0.5 no, to some baryon density n^. The relaxation in this 
zone is provided by the modified Urea processes (b). The 
intermediate zone "D" extends from the density no, at 
which the direct Urea process becomes open, to some den- 
sity nn- The relaxation goes much faster there, because it 
is provided by the direct Urea processes (c) while the mod- 
ified Urea processes are always insignificant once the direct 
Urea processes are open. The inner zone "H" extends from 
the density tih, at which hyperonic processes (d) start to 
operate, to the stellar center. In this zone, we have the re- 
laxation channels of two very different types: direct Urea 
(c) and nonlepton hyperonic processes (d), with th <td. 
The low-mass neutron stars contain only zone M. More 
massive stars may contain zones D and H. 

In a pulsating star, the baryon number density rtb 
varies with pulsation period 2tt/uj ~ 10~ 3 s, where u> ~ 
10 4 s _1 is a typical pulsation frequency. Even if the unper- 
turbed (non-pulsating) stellar state was in the full equi- 
librium (which we assume), the pulsating matter may be 
out of the equilibrium with respect to channels (b)-(d). 
The adiabatic index for perturbations of the matter has, 
therefore, to be defined as 



P \dn h 



P/c 2 



P 



dP 



(9) 



where p is the mass density and the subscript C means 
that the derivative should be constrained by the conditions 
in various channels (see below). In the full equilibrium, we 
have T\ = Teq, where Teq is the adiabatic index that de- 
termines the stiffness of the equation of state (EOS) in 
a non-pulsating star (i.e., determines the basic neutron- 
star structure in equilibrium). Figure 1 shows the density 
dependence o f Teq for a relativistic mean-field EOS of 
Glendenning ( |l985| ; his case 3) of npe^AE - matter (sim- 
ilar curve of Teq was calculated, for instance, by Balberg 
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Fig. 1. Adiabatic index Ti versus baryon number density 
(lower horizontal scale) or mass density (upper horizon- 
tal scale) in zones M, D, and H (indicated by arrows). 
The curves show Leq at full equilibrium and Tfr at fully 
frozen composition of matter. Vertical dotted lines: the 
thresholds of appearance of /x, A, and S~. 



& Gal 1997). For this EOS, fi, A, and S~ particles ap- 
pear at nb = 0.110, 0.310, and 0.319 fm -3 , respectively, 
and the main direct Urea process, Eq. ([|), with electrons 
is allowed at no = 0.227 fm~ 3 . The adiabatic index Ti 
in question may differ from Teq because of the lack of 
equilibrium in channels (b)-(d) in pulsating matter. 

Exact calculation of Ti is complicated. Thus we re- 
strict ourselves to several important limiting cases. 

Generally, the pressure depends not only on nb but also 
on the fractions Xj = ni/n^ of various particle species j. 
Thus, Eq. (|) can be written as 



(10) 



where Tfr = (dlnP/dlnn^)^^ is the adiabatic index 
calculated for the "frozen" composition of matter (at Xj = 
constant). The frozen-composition index Tfr is also de- 
picted in Fig. 1 for the same model of dense matter. 
Freezing the particle fractions makes an EOS stiffer and 
increases Ti. The real value of Ti lies evidently between 
Teq and IVr. Some examples are shown in Fig. 2. 

Two obvious conditions of varying the parameters of 
matter while calculating the derivatives in Eq. ([h]) are: 
conservation of baryon number and electric neutrality. 
This reduces the number of independent fractions (degrees 
of freedom) by two. Other conditions are described below. 

3.1. Subthermal pulsations in nonsuperfluid matter 

Let us consider nonsuperfluid matter and the subthermal 
pulsation regime, Sfij ^ T, for all j. In zones M and D, the 
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Fig. 2. Same as in Fig. |l| but for different equilibration 
regimes. Thick solid line: I\ for subthermal perturbations 
with Urea channels (b) and (c) frozen, but hyperonic chan- 
nels (d) fully relaxed. Thick dash-and-dotted line in zones 
D and H: I\ = T[ for suprathermal perturbations with 
all the channels, but modified Urea, fully relaxed. Thin 
dashed curves: Teq and Tfr as in Fig. [j]. 

relaxation times tm and td are much larger than the pul- 
sation periods 2it/uj. Accordingly, Ti = Tfr is a good ap- 
proximation, i.e., Ti is noticeably higher than Teq (Figs. 1 
and 2). 

In zone H, the situation is more complicated. If the star 
is sufficiently cold, T <C 10 9 K, we have tr 3> 2tt/lu, and 
again 1^ = Tfr. However, in a hotter star, we may have 
the opposite condition, r H <C 2tt/uj, which will lead to the 
equilibrium in hyperonic channels (d); then Ti has to be 
calculated assuming this equilibrium. The latter implies 
two conditions: 



"n + "A=A* P + Ms, 



XT. + XA 



(11) 



in pulsating matter (the upperscript (0) refers to equilib- 
rium values). In combination with conservation of baryon 
number and electric neutrality in pulsating matter, condi- 
tions (O) allow us to calculate the partial derivatives in 
Eq. (|l0|) and determine I\ (plotted by the thick solid line 
in Fig. 2) . This adiabatic index is only slightly lower than 
TpR, for the adopted EOS (cf. the thick solid line and the 
upper dashed line in Fig. 2). 

3.2. Suprathermal pulsations in nonsuperfluid matter 

Let us outline the case of larger pulsation amplitudes, 
6/j,j 3> T, which however are assumed to be relatively 
small, 6/ij <C fij, so that the linear pulsation theory 
is valid. Under such conditions, the deviations of the 
chemical potentials from the fully equilibrium values will 
strongly reduce the relaxation times r. As follows from 
calculations of the rates of nonequilibrium processes (e.g., 
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Reisenegger 1995, Yakovlev et al. 2001), one can estimate 
the relaxation time at Sfij 3> T by replacing 



T —*T = \JT 2 + (<5m)7(2tt) 2 (12) 

in Eqs. (^), where Sfi is a typical pulsation amplitude 
of chemical potentials. This replacement correctly repro- 
duces the regime of subthermal pulsations (Sfi T, 
T sa T). In the suprathermal regime, Sfi 3> T, the relax- 
ation times no longer depend on T but depend on the pul- 
sation amplitude Sfi instead. Let us stress that the replace- 
ment is approximate; detailed calculations are required to 
obtain quantitatively accurate results in the suprather- 
mal regime. Such calculations are beyond the scope of the 
present paper. 

To illustrate the suprathermal effects let us take, for 
instance, T = 100 keV - 10 9 K, /z ~ 500 MeV, and Sfi ~ 
0.02 fi ~ 10 MeV. Then we have Sfi/T ~ 100 and the 
relaxation times tm ~ 0.1 s, td ~ 10 -4 s, and th ~ 3 x 
10 -6 s. Thus, tm remains much larger than 2tt/uj, and we 
have Ti = Tfr, as before, in zone M. However, td.h <S 
2~k jio, and we have rapid relaxation in all channels in zones 
D and H, where we can put T± = Teq. The adiabatic index 
for this regime is plotted by the dot-and-dashed line in 
Fig. 2 (marked as T[). 

3.3. Subthermal pulsations in superfluid matter 

Now we return to subthermal pulsations, Sfi ^ T, but as- 
sume superfluidity of baryons. All baryon species (n, p, A, 
and E~, in our example) can be in a superfluid state. The 
critical temperatures T c of baryon superfluids depend on 
the model of strong interactions and the many-body the- 
ory employed. Various microscopic models predict density 
dependent T c , different for different baryons, rangin g from 
- 10 8 K to - 10 10 K (e.g., Lombardo & Schulze pOOl 



also see Haensel et al. 2001 , 2002 for references). A strong 
superfluidity (T c 3> T) will drastically (exponentially) re- 
duce the rates of the processes involving baryons (due to 
the appearance of superfluid gaps in baryon energy spec- 
tra; see, e.g., Yakovlev et al. 2001). This will increase the 
relaxation times (||) in various channels (b)-(d). In this 
case we can easily have r 3> 2tt/lu, and the relaxation will 
be frozen over pulsation periods. Thus, putting Ti = Tfr 
will be a good approximation. 

3.4. Suprathermal pulsations in superfluid matter 

This case is complicated. Note however that the pulsations 
with amplitudes Sfi higher than the superfluid gaps, A 
1 MeV, seem to be independent of superfluid properties 
of dense matter. Thus, the relaxation times of pulsations 
with such amplitudes can be estimated using the results 



of Sect. 3.2 



4. Damping of stellar pulsations by bulk viscosity 

Let us turn to the bulk viscosity of dense matter and asso- 
ciated damping of neutron-star pulsations. For the compo- 



sition of dense matter we are interested in, these problems 
have been studied in a number of papers starting from the 
pioneering papers by Langer & Cameron (1969) and Jones 
(1971). The current state of the problem and references to 
other papers can be found in Jones (2001), Haensel et al. 
( p002| ), and Lindblom & Owen fl2002| ). 



4.1. Nonsuperfluid matter 

Simple estimates of the bulk viscosity £ in the familiar sub- 
thermal regime are given, for instance, by Haensel et al. 
( 2001 . [2002 ). The generalization to the case of suprather- 
mal pulsations is straightforward. ^Generally, one can in- 
troduce an effective bulk viscosity £ which determines ac- 
tual dissipation rate of pulsation energy per unit volume. 
The estimates of £ can be obtained from the estimates of 
C by replacing there T — > T. In this way we come to the 
order-of-magnitude estimates of the effective bulk viscos- 
ity in zones M, D, and H: 



Cm - 5 x 10 18 w 4 - 2 T 9 e g cm" 1 s" 1 , (13) 

Cd ~ 5 x 10 24 ^ 2 T 9 4 g cm" 1 s-\ (14) 

Ch ~ 10 29 W4 - 2 T 9 2 gem" 1 s-\ (15) 

where u>4 is the pulsation frequency in units of 10 4 s~ 4 . 

At the next step, following Haensel et al. (2002), we 
estimate a typical viscous dissipation time of stellar pul- 
sations, idiss ~ pR 2 1 'C, with R ~ 10 km being the radius 
of the neutron-star core. The timescale idiss is determined 
either by dissipation in zone M (idiss = iivi, if zones D 
and H are absent), or dissipation in zone D (idiss = tr>, 
if zone H is absent), or dissipation in zone H (idiss = iH, 
if a representative zone H is available). Then we obtain 
drastically different dissipation times, 



hi 



lOoj 



yr, t u 



100 uj\ 



1 Q 



ms 



(16) 



for the different cases. 

For the subthermal pulsations, we have T — T and re- 
produce thus the estimates given by Haensel et al. (2002). 
For the suprathermal pulsations, Sfi S> T, the effective 
bulk viscosity and dissipation time (just as the relaxation 
times in Sect. 3.2) stop to depend on temperature T but 
depend on Sfi instead, producing strongly enhanced non- 
linear (in the pulsation amplitude, Sfi) viscous pulsation 
damping. The damping will completely destroy the pulsa- 
tion mode if the dissipation time idiss becomes comparable 
to the pulsation period 2tt/uj. Thus, setting idiss = 2tt/lo 
we can find the maximum pulsation amplitude <5/^ max 
which a pulsating star can sustain. For the pulsations in 
regimes M, D, and H, we obtain: 



Sfi 



M 

max 



50wl /2 MeV, 



(17) 
(18) 

, 4 xv.cv. (19) 

Accordingly, the viscous damping cannot allow the star to 
pulsate with too large amplitude, and the maximum am- 
plitude does depend on the damping regime. It is natural 



10w 4 /4 MeV, 
2u 3 J 2 MeV. 
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to assume that the strong damping with the maximum 
amplitude will rapidly reduce the pulsation amplitudes to 
a lower level with a smaller damping rate. 

4.2. Superfluid matter 

The effects of superfluidity are similar to those described 
in Sects. 3.3 and 3.4. For subthermal pulsations, these 
effects were considered, for instance, by Haensel et al. 



( |200fj| , |2001| , |2002| ). A strong superfluidity may drastically 
(by many orders of magnitude) reduce the bulk viscos- 
ity and the dissipation rate. The effects of superfluidity 
on the bulk viscosity and the damping of pulsations in 
the suprathermal regime with 5/i i> A are expected to be 
small, i.e., the results of Sect. 4.1 can be used. 



5. Summary 

We have outlined the most important factors which affect 
the adiabatic index Ti for density perturbations and the 
viscous damping of pulsations in the neutron-star cores. 
First, non-leptonic strangeness-changing processes with 
hyperons may be rapid enough to establish partial hy- 
peronic equilibrium over pulsation periods and decrease 
Ti slightly below the "frozen" adiabatic index, Tfr, in 
hyperonic matter. In addition, they induce a rapid vis- 
cous damping even in the subthermal pulsation regime 
(5[ij ^ T). Second, suprathermal but still linear pertur- 
bations (T <C Sfij <C fij) may produce a very rapid equili- 
bration in various relaxation channels which reduces Ti to 
the fully equilibrium value Teq , enhances the bulk viscos- 
ity and damps the pulsations quickly to lower amplitudes. 
Third, superfluidity of baryons in the subthermal regime 
increases relaxation times in various channels, brings T± 
closer to Tfr, and reduces the viscous dissipation. 

Therefore, a proper calculation of neutron star pulsa- 
tions and their dynamical evolution represents a compli- 
cated problem. Much work is required to study this prob- 
lem in full detail. Our order-of-magnitude estimates have 
to be replaced with accurate numerical solutions of the 
equations of stellar pulsations taking into account proper 
boundary conditions and joint effect of various factors in 
all neutron-star layers, from the surface to the center. Our 
assumption of one typical pulsation amplitude of chemi- 
cal potentials, (5/i, in Eq. ( fl^ ) is an oversimplification. In 
reality, one has to deal with the number of density depen- 
dent amplitudes 5/j.j for different particle species j. The 
relaxation in different equilibration channels and the vis- 
cous damping of pulsations can be strongly nonuniform; 
one cannot exclude the existence of thin layers in the neu- 
tron star cores, where the damping is exceptionally strong. 
For instance, they may be the layers where new hyperons 
appear, with sufficiently small chemical potentials Hj just 
beyond their appearance threshold. Neutron star pulsa- 
tions in these layers may be suprathermal, producing en- 
hanced damping. The problem of suprathermal pulsation 
regime, without and with superfluidity of baryons, is of 



special importance. This regime will be accompanied by 
huge energy release which will heat the star. 
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